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ABSTRACT 

We present a numerical model for the extended steep-spectrum radio sources and 
the elongated X-ray structure in A3667 based on new 3-dimensional MHD/N-body 
simulations. The X-ray and optical analyses of A3667 indicate that it has undergone a 
recent subcluster merger event. We believe that the Mpc-scale radio sources identified in 



A3667 are also a consequence of the merger. Our previous numerical simulations show 
that mergers often produce large-scale shocks and turbulence capable of both magnetic 
field amplification and in-situ reacceleration of relativistic particles. Our model suggests 
that these radio structures, separated by ~2.6/ifoo Mpc, are in fact causally linked via 
a slightly off-axis merger that occurred nearly in the plane of the sky approximately 1 
Gyr ago with a subcluster having a total mass equal to ~20% of the primary cluster. 



Subject headings: magnetohydrodynamics - methods: numerical- galaxies: intergalac- 
tic medium - galaxies: clusters: individual (Abell 3667) 



1. INTRODUCTION 

Abell 3667 (see Fig. HD at z=0.055 (Sodre et 
al. 1992) is a massive (> 10^^ M©; Sodre et 
al. 1992; Knopp et al. 1996), X-ray luminous 
(~10'^^/i~^ erg s~^; Knopp et al. 1996) cluster 
of galaxies which exhibits a variety of the ob- 
servational signatures of a recent cluster merger 
(Roettiger et al. 1996). Here, we present a self- 
consistent numerical MHD/N-body model sug- 
gesting that the most unique features of A3667, 
the Mpc-scale radio sources (Rottgering et al. 1997), 
are similarly a signature of a recent merger. 

Characterized variously as a "halo" or a "relic" , 
we believe that the extended radio emission in 
A3667 represents a stage in the evolution of a 
radio halo. Radio halos are generally described 
as large (>0.5 Mpc), amorphous steep-spectrum 
structures which do not appear to be associated 
with any particular galaxy and, as such, are be- 
lieved to be intrinsic to the cluster itself (Hanisch 
1980). By contrast, a radio relic was once asso- 
ciated with a particular AGN that has turned- 
off and possibly drifted away leaving the radio 
source to diffuse and fade away. It is possible 
that radio relics supply the seed relativistic par- 
ticles necessary to produce radio halos, but the 
connection is unclear at this time and we do not 
address this issue here. 

Theoretically, the difficulty in understanding 
radio halos is related to their large spatial ex- 
tent. If the relativistic electrons responsible for 
the radio synchrotron emission originated from 
a single source (e.g., AGN), they would be re- 
quired to diffuse over Mpc-scales during their rel- 
atively short radiative lifetimes (of order 10® yrs) 
requiring velocities well in excess of the limiting 
Alfven speed. (Holman, lonson & Scott (1979) 
present an alternate view.) This has led sev- 
eral researchers to suggest that radio halos re- 
quire a means by which relativistic particles can 
be accelerated in-situ, such as via galactic wakes 
(Roland 1981) or shocks and MHD turbulence 
(e.g., Eilek & Henriksen 1984) generated during 



a cluster merger (DeYoung 1992; Tribble 1993b). 
Whereas galactic wakes, appear to be insufficient 
to power the radio halo (Goldman & Repheali 
1991; DeYoung 1992), cluster mergers, having a 
kinetic energy comparable to the total thermal 
energy of the cluster, are more than sufficient to 
power the radio source (Tribble 1993b; Burns et 
al. 1994; Burns 1998). Mergers may also over- 
come part of the electron transport problem by 
supplying bulk flows (>1500 km s~^ ; Roettiger 
et al. 1997) capable of transporting the relativis- 
tic electrons over large distances during their ra- 
diative lifetime (Rottgering et al. 1994). 

Observationally, mergers have been linked to 
radio halos by both a correlation with substruc- 
ture (X-ray and optical) and by an anticorrela- 
tion with cooling flows (Edge, Stewart & Fabian 
1992). For example, Coma, A2255 and A2256 
all contain radio halos of one type or another 
and exhibit substructure that is indicative of re- 
cent dynamical evolution. It has also been sug- 
gested that cooling flows (e.g., see Fabian 1994 
for a review of cooling flows) will be disrupted 
by massive mergers (McGlynn & Fabian 1984; 
Edge et al. 1992; Burns et al. 1997). Gomez 
et al. (1999), using numerical simulations, have 
shown that mergers can disrupt cooling flows un- 
der various initial conditions. None of the above 
mentioned clusters contains a cooling flow, nor 
for that matter does A3667. 

The purpose of this paper is to demonstrate 
using fully 3-dimensional, numerical MHD/N- 
body simulations that the X-ray, optical and ra- 
dio observations of A3667 are consistent with a 
single merger model in which shocks formed dur- 
ing the merger provide sites for diffusive shock 
acceleration of cosmic ray (CR) electrons (see 
Longair (1994), and references therein, for a re- 
view of particle acceleration in shocks). We do 
not attempt to model the CR acceleration, dif- 
fusion, advection, and subsequent "aging" self- 
consistently with the MHD. When done in de- 
tail, this is a difficult problem and can be com- 
putationally expensive. Recently, Jun & Jones 



(1998) and Jones, Ryu & Engler (1998) presented 
a promising approach to this problem (apphed to 
supernova remnants) which we will look to incor- 
porate into future work. At this time, we present 
a model in which we introduce a population of 
CR electrons based on the shock characteristics 
during an epoch believed to be representative of 
A3667's current dynamical state, as inferred from 
the X-ray and optical data. The CR diffusion, 
advection, and aging are modeled by assuming 
that the MHD properties (i.e., shock strength, 
structure and velocity, bulk flows etc.) do not 
evolve significantly during the radiative lifetime 
of the particles (~10* years) and that particle 
reacceleration behind the shock is unimportant. 
The modeling presented here is similar to that of 
Tribble (1991, 1993a, 1994) but with two impor- 
tant exceptions. His models were purely static 
while ours are based on the self-consistent evo- 
lution of shocks and magnetic fields within the 
merger model, and we are attempting to specifi- 
cally address the characteristics of the extended 
radio emission in A3667. Although others have 
suggested that A3667 has recently undergone a 
merger (Knopp et al. 1996), and still others have 
suggested a possible shock origin for the extended 
radio emission (Rottgering et al. 1997; Ensslin 
et al. 1998), this is the first study that presents a 
quantitative numerical model which attempts to 
simultaneously explain the radio, optical and X- 
ray morphology of this cluster within the context 
of a single merger. 

Section ^ is a review of recent observations 
and evidence for substructure in A3667. In g^, 
we discuss our numerical method. Our initial 
conditions are presented in §Q. Section ^ sum- 
marizes our model constraints. We discuss the 
radio halo model and its limitations in §^. Fi- 
nally, we summarize our results in ^ The model 
is scaled assuming Ho=70 km s~^ Mpc~^. The 
observational data (see §|2|) is parameterized by 
h=Ho/100. 

2. A REVIEW OF THE OBSERVATIONS 



2.1. The X-ray Data 

Knopp et al. (1996) provide a detailed analy- 
sis of the ROSAT X-ray data which we summa- 
rize here. The X-ray surface brightness (XSB; 
contours, Fig. Q) is seen to be elongated from the 
southeast (SE) to the northwest (NW). The ellip- 
ticity is greatest in the cluster core and decreases 
with increasing radius. There is an isophotal 
twisting of ~15° from the inner to outer cluster. 
The XSB exhibits an extension to the NW of the 
core while it is seen to be flattened to the SE. 
This results in a significant centroid shift (Mohr 
et al. 1995). A3667 does not appear to be isother- 
mal. Knopp et al. find a mean temperature of 
6.3±Q;g keV. However, the core is found to be 
hotter (9±f:[] keV) than the mean at the 82% 
confidence level. Markevitch et al. (1998), using 
ASCA data, find a mean temperature of 7.0±0.6 
keV. Using the ROSAT-derwed XSB, Knopp et 
al. find the total X-ray luminosity within 0.80/i~^ 
Mpc to be L^=6.74±0.09x lO'^^/i-^ erg s'^. They 
estimate the gas mass within 0.80/i~^ Mpc to be 
8.9±0.9xl0^2 /j-5/2 jyj^ ^j^-jg ^Yie total dynam- 
ical mass within the same radius (assuming hy- 
drostatic equilibrium, which is questionable in a 
recent merger) is 3.5ibQ4xlO^^ h~^ M©, giving a 
relatively low baryon fraction of 2.5/i~^'^%. The 
total mass within 3.0 h~^ Mpc is estimated to be 
1.4 ±0.2x 10^^ /i"^ M0. 

2.2. The Optical Data 

Sodre et al. (1992) provide 203 galaxy posi- 
tions (128 redshifts) in the vicinity of A3667. The 
mean cluster redshift is 16557ibll7 km s^^, and 
the velocity dispersion is 1223ib72 km s~^. Fadda 
et al. (1996), using a somewhat larger galaxy 
sample and an analysis of substructure, find a 
global velocity dispersion of 971±4y although the 
velocity dispersion in the core appears to be 
much higher. Using the global velocity disper- 
sion, and the ASCA mean temperature, one cal- 
culates /3spec=0.84. This value is somewhat dis- 
crepant from the value of /3/it =0.535 (Knopp et 



al. 1996). Numerical simulations have linked the 
so-called "/3-discrepancy" (Sarazin 1988) to re- 
cent merger activity (Navarro, Frenk & White 
1995; Roettiger et al. 1996). Sodre et al. fur- 
ther note two spatially distinct clumps within the 
galaxy distribution. One clump is centered on 
the brightest D galaxy and the XSB peak while 
the other is associated with the second bright- 
est D galaxy located ~16.6' to the NW (but 
still within the observed XSB distribution). Fig- 
ure |2| shows the positions of galaxies having red- 
shifts supplied by Sodre et al. and the results of 
the (5-test for substructure (Dressier & Shectman 
1988). The circles, marking individual galaxy po- 
sitions, are scaled in radius according to their re- 
spective 5 parameter. The value of 6 indicates 
the degree to which a galaxy and it's 10 near- 
est neighbors deviate from the global velocity 
properties. In Figure |2[ we have filled the cir- 
cles corresponding to galaxies having the lowest 
6 parameters, thereby accentuating the location 
of the two groups noted by Sodre et al. . The 
degree of substructure demonstrated by the 6- 
test gives some indication as to why the observed 
range of velocity dispersions, and therefore clus- 
ter virial masses, is so large. Sodre et al. (1992) 
estimate a total mass of 2.6x10^^ h~^ Mq for the 
entire cluster while Biviano et al. (1993) find 
1.6x10^5/1"^ Mq within 0.75/1"! j^p^. Both val- 
ues are larger than the X-ray derived values (see 



2.3. The Radio Data 

There are many radio point sources located 
within the spatial extent of A3667, but for our 
purposes, we are concerned only with the large 
extended radio sources located to the SE and 
NW of the X-ray core (See Fig. |^, grayscale). 
The uniqueness of the extended radio sources in 
A3667 have been known for many years (Schillizzi 
& McAdam 1975; Goss et al. 1982; Jones & 
McAdam (1992)). Here, we summarize the most 
detailed analysis to date, Rottgering et al. (1997). 
Approximately 1.7h~^ Mpc NW of the cluster's 



X-ray core (and beyond the extent of the XSB), 
is an extremely large {r^lAh~^ Mpc) arc-shaped 
source that does not appear to be associated with 
any particular galaxy in the cluster. Along the 
NW rim of this source, there is a sharp edge 
to the synchrotron distribution. The source be- 
comes more diffuse to the SE. The spectral in- 
dex appears to be relatively flat (a=0.5; S^ oc 
I/"") along the NW edge and exhibits a gradi- 
ent toward the SE steepening to a=1.5. Overall, 
a ~1.1. 

To the SE of the X-ray core there are two 
large, linear sources situated perpendicular to 
the X-ray major axis (parallel to each other). 
The largest and furthest SE extends ~0.7/i~! 
Mpc and is ~0.7/i~! Mpc from the core. The 
source closest to the core has been identified 
by Hunstead and Wieringa (1998) as a narrow- 
angle tailed source (NAT) associated with galaxy 
157 from Sodre et al. (1992) and is therefore 
not of interest to this study. The NW and SE 
halo sources are separated by nearly 2.6h~^ Mpc. 
Both sources are located near or beyond the ex- 
tent of i^O^^T-observed X-ray emission. 

3. NUMERICAL METHOD 

The numerical method is the same as that 
used in Roettiger, Stone & Burns (1998b). The 
ICM and magnetic fields are evolved using ZEUS 
(Stone & Norman 1992a, b), an Eulerian, finite- 
difference code which solves, self-consistently, the 
equations of ideal MHD. The numerical evolution 
of the magnetic field components is performed by 
the constrained transport (CT) algorithm (Evans 
& Hawley 1988) which guarantees preservation of 
the divergence- free constraint at all times. The 
method of characteristics (MOC) is used for com- 
puting the electromotive force (Hawley & Stone 
1995). An extensive series of MHD test problems 
have demonstrated that the MOC-CT method 
provides for the accurate evolution of all modes of 
MHD wave families (Stone et al. 1992). We em- 
ploy outflow boundary conditions on the MHD. 



The collisionless dark matter is evolved using 
an N-body code based on a standard particle- 
mesh algorithm (PM, Hockney & Eastwood 1988). 
The particles and gas are evolved on the same 
grid using the same time step. The time step 
is determined by applying the Courant condi- 
tion simultaneously to both the dark matter and 
the magnetohydrodynamics. The only interac- 
tion between the collisionless particles and the 
gas is gravitational. Since we are modeling an 
isolated region, the boundary conditions for Pois- 
son's equation are determined by a multipole ex- 
pansion (Jackson 1975) of the total mass distri- 
bution (dark matter and gas) contained within 
the computational grid. Particles that leave the 
grid are lost to the simulation. Typically, less 
than a few percent of the particles leave the grid. 

The hybrid ZEUS/PM code was parallelized 
using the Message Passing Interface (MPI; Gropp, 
Lusk & Skjellum 1994). These simulations were 
run on the Cray T3E in the Earth and Space 
Data Computing Division of the NASA Goddard 
Space Flight Center. 

The simulation whose results are presented 
here had an effective resolution scaled to 18 kpc, 
or more significantly, there are 16 zones per pri- 
mary cluster core radius. The computation was 
performed on a grid having dimensions 512 x 
256^. The grid is uniform from zone 100 to 412 
along the merger axis and from zone 90 to 166 
perpendicular to the merger axis. Outside of this 
central region, we gradually increase the zone size 
by ~3% from one zone to the next. 

4. INITIAL CONDITIONS 
4.1. Dark Matter and Gas 

Our initial conditions are similar to those used 
in our previous studies (e.g., Roettiger, Burns 
& Loken 1993; Roettiger, Loken & Burns 1997; 
Roettiger, Stone & Mushotzky 1997, 1998a). We 
begin with two clusters whose gas distributions 
are consistent with observations of relaxed sys- 
tems. Our cluster dark matter distributions may 



be less consistent with observations in that they 
do not have the central density cusp implied by 
recent strong lensing observations and numeri- 
cal simulations of large-scale structure formation 
(Navarro, Frenk, & White 1997). We do not be- 
lieve that this difference will significantly alter 
the merger dynamics. The mass distributions in 
our simulations are based on the lowered isother- 
mal King model described in Binney & Tremaine 
(1987). The lowered isothermal King model is a 
family of mass distributions characterized by the 
quantity V'/c^ which essentially defines the con- 
centration of matter. As ip/cr'^ increases, the core 
radius {re) decreases with respect to the tidal ra- 
dius (rt). We have chosen a model with ip/a'^ =12 
in which we have truncated the density distribu- 
tion at 15rc. Near the half-mass radius, the to- 
tal mass density follows a power-law distribution, 
p r^ r~", where a ^ 2.6. This model is consis- 
tent with mass distributions produced by cosmo- 
logical N-body simulations which show a ^ 2.4 
in a high density universe (fi=l) and --^2.9 in 
a low density universe (r2=0.2) (Crone, Evrard, 
& Richstone 1994). Initially the gas distribu- 
tion (ICM) is in hydrostatic equilibrium within 
the gravitational potential defined by both the 
gas and dark matter components. The gas dis- 
tribution is isothermal within the central 6rc- 
At larger radii, the temperature drops gradually. 
The density peaks are initially separated by 6.5 
Mpc. Each cluster was given a small initial ve- 
locity in order to speed up the merger process 
and thus conserve computational resources. The 
final impact velocity (~2750 km s~^ ) is not af- 
fected significantly by the initial velocity which 
has components of ~275 km s~^ parallel to the 
line of centers and -^30 km s^^ perpendicular to 
the line of centers, resulting in a slightly off-axis 
merger. The physical scaling of the pre-merger 
clusters can be found in Table 1. 

4.2. Magnetic Field 

The only galaxy clusters certain to have large- 
scale or intrinsic magnetic fields are those with 



radio halos, and these sources are rare (Hanisch 
1982). Of those containing halos, Coma has been 
studied in the most detail (e.g., Jaffe 1977; Kim 
et al. 1990 among others). Deiss et al. (1997) 
show that the synchrotron halo in Coma, after 
subtracting point sources, traces the X-ray sur- 
face brightness distribution indicating that the 
magnetic pressure gradient is similar to that of 
the thermal pressure. Of course Coma is only 
one cluster and may not be representative of ha- 
los in general. In addition to the spatial distri- 
bution of the magnetic pressure, we must also be 
concerned with the spectral power distribution. 
Cluster magnetic fields are believed to exist as 
flux ropes that are tangled on a variety of phys- 
ical scales (e.g., Ruzmaikin et al. 1989). Ob- 
servations of polarization in cluster radio sources 
indicate tangling occurs on scales ranging from 
a kiloparsec (e.g., Feretti & Giovannini 1997) to 
tens or even hundreds of kiloparsecs (Kim et 
al. 1990; Kim, Kronberg, & Tribble 1991; Feretti 
et al. 1995) . Studies of cluster radio halos using 
a variety of methods (see Miley 1980) indicate 
that typical mean field strengths are of order 1 
//G (e.g.. Coma, Kim et al. 1990; Ensslin & Bier- 
mann 1998; A2255, Burns et al. 1995; A2256, 
Rottgering et al. 1994; Bagchi, Pislar & Lima 
Neto (1998)). 

In order to model a random field tangled 
on a variety of scales, we begin by defining a 
magnetic vector potential, A(k)=Aok~ , where 
the amplitudes (Ao) for each Cartesian coordi- 
nate are drawn from a Gaussian distribution. 
The vector potential is then transformed via 
a 3-dimensional Fast Fourier Transform (FFT) 
into physical space where it is scaled spatially 
by the gas density distribution. Assuming a 
uniform spherical collapse (likely a gross over- 
simplification, e.g., Evrard 1990; Bryan et al. 1994; 
among many others) and flux freezing, it can 
be shown that the magnetic field will scale as 
pjas or B'^ r^^'^. Finally, we initialize a tan- 
gled divergence-free magnetic field from A via 
B = V X A. This results in B^ oc \?^^-^^ . Here 



we adopt 5 = 5/3. Unfortunately, neither ob- 
servation or theory are capable of constraining 
our choice of power-law index. We have scaled 
the magnetic pressure such that the mean field 
is a dynamically insignificant 0.3 [iG within 2rc. 
Due to numerical resolution limitations, the min- 
imum resolved magnetic field scale is ~4 zones or 
nearly 80 kpc. 

5. THE MODEL 

Based on our numerical model and limited sur- 
vey of merger parameter space (cluster mass ra- 
tios, impact parameters, relative core densities, 
etc.), we believe that A3667 has undergone a 
merger approximately 1 Gyr ago with a subclus- 
ter having a mass 20% of the primary cluster. 
The merger was slightly off-axis having an im- 
pact parameter of <0.5rc. The merger occurred 
largely in the plane of the sky with the subclus- 
ter moving from the SE to the NW passing on 
the southern side of the more massive or pri- 
mary cluster. Figure ^ shows the synthetic XSB 
(dashed contours) overlaid onto the projected, 
emission-weighted temperature (solid contours) 
at the epoch determined to most closely resemble 
A3667. The sharp edges in the temperature dis- 
tribution indicate the locations of shocks which 
we associate with sites of particle acceleration 
for the purposes of modeling the radio emission 
(See §^). Figure ^3 is a comparison of the sim- 
ulated and ^^C^-observed temperature profiles 
(Markevitch et al. 1998). Figure ^ shows the 
synthetic XSB (solid contours) overlaid on the 
simulated synchrotron emission (grayscale; see 
§g). The -|- symbols indicate the location of dark 
matter /galaxy centroids. The primary cluster is 
largely coincident with the XSB peak while the 
subcluster remnant is located to the NW. The fu- 
ture evolution of this system would be similar to 
that of A754 as described in detail by Roettiger 
et al. (1998a). 

In producing our model of A3667, we have 
made some effort to match as many of A3667's 
observational characteristics (see §||, Figs. || and 



^ as closely as possible, while keeping in mind 
that our primary concern for this work is the 
merger- induced shock structure which, as our 
previous simulations show (e.g., Roettiger et 
al. 1997) are a rather ubiquitous feature of ma- 
jor merger events. Regarding the XSB distribu- 
tion, we are able to reproduce the overall elonga- 
tion which defines the merger axis, the isopho- 
tal twisting which indicates a slightly off-axis 
merger, and the steepness of the XSB distribu- 
tion to the SE relative to the NW. We believe 
that the subcluster first impacted from the SE 
moving to the NW causing a compression of the 
X-ray emitting gas on the SE side. As the sub- 
cluster dark matter remnant exits the primary 
cluster core, it draws out the ICM creating the 
NW extension. With regard to the galaxy distri- 
bution, we are able to reproduce the relative po- 
sitions of the galaxy concentrations (dark matter 
particles in our simulations, +'s in Fig. |^) with 
respect to the XSB and the radio sources. Below, 
we describe radio source modeling in greater de- 
tail below. 

6. DETAILS OF THE RADIO HALO 
MODEL 

6.1. Basic Assumptions 

There are several simplifying assumptions in- 
herent in this model. First, we assume that 
all particle acceleration occurs at the shocks. 
Once the particles leave the region of the shock, 
there is no further acceleration. Second, we as- 
sume a simple diffusive shock model in which the 
strength of the shock (as measured by the density 
compression) determines the power-law slope of 
the injected particle energy spectrum (e.g., Drury 
1983; Longair 1994). The electrons are assumed 
to be accelerated within a single zone. That is, 
their scattering length is small compared to a 
computational zone, and the time required to ac- 
celerate them is small compared to a time step. 
Third, we assume that the basic shock structure 
(strength, shape, velocity, etc.) does not change 



significantly during the radiative lifetime of the 
relativistic electrons. Finally, we assume that 
some small fixed fraction of the thermal particle 
flux is accelerated to its equilibrium distribution 
as it passes through the shock. 

6.2. Method 

Having chosen the epoch most representative 
of A3667 based on the X-ray and optical data, 
we begin modeling the radio emission by identi- 
fying the locations of the shocks. In ZEUS, an 
artificial viscosity is used to smooth shocks over 
4-5 zones. We locate the shock by calculating 
the gas compression ratio (r = P2I P\) at two 
points separated by 5 zones throughout the com- 
putational volume. The compression ratio can 
be related directly to the power-law index (7) 
of the injected relativistic electron energy spec- 
trum, N{E) ex E-'^ such that 7 =^^+1. Once 
identified, each point along the shock is assigned 
a value of 7 based on the local strength of the 
shock. For the most part, the shock strength is 
relatively uniform across the face of the shock, 
however it does weaken somewhat near its edges. 
The shock structure present within the plane of 
the merger can be seen in Fig. |5|. Note that 
there are multiple shocks to the SE of the X-ray 
core each presenting a potential site for particle 
acceleration. 

The next step is to associate each zone behind 
(or upstream from) the shock with the nearest 
point on the shock along the fluid flow lines. We 
then assume that the particles in a given zone 
behind the shock are derived from that location 
on the shock and that their distance (d) from 
the shock is directly proportional to the veloc- 
ity of the shock relative to the bulk fluid flow 
(Vg), and the time since acceleration (to) such 
that d = nVsta- The scaling factor k is used 
to parameterize the diffusion/advection rate. If 
the electrons diffuse freely in straight-line mo- 
tion, unaffected by the magnetic field or turbu- 
lence, K=l. This would represent the limiting 
case in which the magnetic field is weak and/or 



aligned parallel to the flow and accelerated par- 
ticles are simply deposited in the wake of the 
shock. Of course, the situation is almost cer- 
tainly more complicated. At the other extreme 
(strong, tangle field), Tao (1995) suggests that 
the maximum inhibition of electron diffusion by 
magnetic fields is of order a factor of ten. That 
is, the effective path length of an electron trav- 
eling through a tangled magnetic field environ- 
ment will be at most ten times longer than the 
straight-line path, (i.e., k=0.1). Here, we adopt 
the weak field/high diffusion limit where k ~1.0. 

In order to quantify the rate at which particles 
diffuse/advect from the shock in the absence of a 
magnetic field (i.e., Vs), we have incorporated a 
passive-scalar field into ZEUS. The passive-scalar 
is a dynamically insignificant quantity (in that 
sense, it is likely similar to the radio plasma) 
which is advected with the thermal gas density. 
In this implementation, we inject the scalar field 
at the locations of shocks and simply allow it to 
evolve. As expected, the passive-scalar field (or 
relativistic electrons) is seen to trail the shocks 
as they propagate through the ambient medium. 
This is consistent with observations of shocks 
within the solar system (Kennel et al. 1986). 
By evolving the scalar field over a limited time 
interval near the A3667 epoch, we can estimate 
the effective shock velocity for each point along 
the shock front. In this manner, we find effec- 
tive shock velocities of 700-1000 km s~^ while 
the SE shocks have an effective velocity of ~250 
km s~^. In fact, the shock nearest the X-ray core 
is effectively stationary. This analysis shows lit- 
tle evidence of mixing between the relativistic 
and thermal plasmas behind the NW shock. Of 
course this result could be strongly resolution de- 
pendent. It may also depend on the initial condi- 
tions. The shocks in this model, particularly to 
the NW, propagate into a largely non-turbulent, 
homogeneous medium which is certainly not the 
case. 



6.3. Synchrotron Aging, Emissivity, and 
Spectral Index Distribution 

We use the formalism of Myers & Spangler 
(1985) to age our synchrotron spectrum. Their 
work is based on the models of Kardashev (1962) 
and Pacholczyk (1970) and Jaffe & Perola (1973). 
Here, we specifically employ the modification in- 
troduced by Jaffe & Perola which assumes that 
the electrons become isotropized on timescales 
short compared to their radiative lifetimes and 
therefore can be represented by a time-averaged 
distribution of pitch angles. Adiabatic loses are 
not included in this model. Based on our analysis 
of the simulation dynamics and the Rottgering 
et al. (1997) analysis showing the radio source to 
be greatly underpressurized with respect to the 
thermal plasma, we believe that expansion losses 
can be neglected. 

Given the initial electron energy spectrum 
power-law index (7), the magnetic field strength 
(-B), and the time since acceleration, we can de- 
termine the spectral index and relative inten- 
sity of emission between 1.4 GHz and 4.9 GHz. 
The emissivity within a given zone is then sim- 
ply, eu oc {Bsin{6))'^^^v~°', where 9 is the an- 
gle between the magnetic field and the observer's 
line-of-sight. The Myers & Spangler model does 
not account for varying magnetic field strengths. 
Therefore, the field used to age the particles 
is the mean field experienced during the elec- 
tron's lifetime (i.e., the mean magnetic field along 
the line connecting the particle's current location 
and it's presumed point of origin on the shock). 
The magnetic field used to calculate the emissiv- 
ity within a given zone is simply the field within 
that zone. Having calculated the synchrotron 
emissivity within each volume, we perform a line- 
of-sight integration through the emitting volume 
to produce the surface brightness image (Fig. 
|a). 

The rate at which the particles age (or equiv- 
alently, their synchrotron spectrum steepens) is 
determined by the strength of the magnetic field 



traversed by the particles. Although the sim- 
ulation was run with a specified field strength, 
since it was dynamically insignificant, we can 
treat it as a free parameter (within a reason- 
able range of values) for purposes of modeling 
the radio spectral index distribution. Examining 
Fig. 3 of Rottgering et al. (1997), we see that 
the spectral index of the NW source steepens to 
~1.5 at about 0.7 Mpc (Ho=70 km s"^ Mpc"^) 
upstream from the shock. Assuming the Jaffe- 
Perola model and high diffusion (k = 1; see § |6.2D , 
we can reproduce the distribution of spectral in- 
dex if the mean magnetic field within the region 
of the source is ~0.6 fj,G. If we restrict the rate of 
diffusion (i.e., k < 1), the magnetic field can be 
decreased accordingly because particles will take 
longer to diffuse/advect a distance of 0.7 Mpc 
and will have aged spectrally in the mean time. 
It is not possible to use this method to constrain 
the field strength in the SE radio sources. We 
do not have spectral index data, and the shocks 
have not moved significantly during the particle's 
lifetimes. Therefore, we scale the global field ac- 
cording to our analysis of the NW source. As- 
suming we had the correct initial magnetic pres- 
sure distribution (see ^.2|), the mean field in 
the region of the SE radio source is ~1.5 //G. 
The greater field strength results from the SE 
source being closer to the X-ray core than the 
NW source (remember B~ r~^'^, initially), and 
from the greater compression (and amplification) 
of the field on the SE side of the cluster as a result 
of the merger. Again, if we have overestimated 
the electron diffusion/advection rate, as is likely, 
these field estimates will decrease. 

6.4. Limitations of the Model 

Of course the model has its limitations which 
are owed primarily to the large number of free 
and not necessarily independent parameters such 
as relative cluster masses, the shape of their ini- 
tial mass distribution, impact parameter, initial 
velocities, magnetic field distribution, projection 
with respect to the observer, etc. It is compu- 



tationally expensive and physically impossible to 
explore all of parameter space. Consequently, we 
present the best fit to the current data given our 
limited survey of the many free parameters. Even 
so, the largest discrepancies between our model 
and the observations appear to be the strength 
of the NW shock and the location of the most 
northwesterly of the SE shocks. The NW shock 
in our model, based on the observed spectral in- 
dex distribution, appears to be too weak. If the 
spectral index along the NW rim of the NW ra- 
dio source (a=0.5) is correct, we need a shock 
with a density compression near 4 (i.e., a strong 
shock). At the chosen epoch, it is typically less 
than 3 which produces a somewhat steeper spec- 
trum (q -^0.95). This may imply that we need to 
consider an earlier epoch since the shock weakens 
as it expands outward. The shock could also be 
strengthened if the impact velocity of the sub- 
cluster where greater (i.e., more massive clusters 
or a greater initial velocity), or if the exiting 
subcluster remnant were interacting with higher 
density gas at large radii or with material falling 
into the cluster from a large-scale filament. This 
highlights another limitation of the model. These 
numerical simulations are non-cosmological. We 
do not attempt to model the external cluster en- 
vironment such as filaments or spherical accre- 
tion. The structure of galaxy clusters at large 
radii (>1.5/i~^ Mpc) is uncertain from an obser- 
vational standpoint, and, consequently, our ini- 
tial conditions become potentially less accurate 
at larger radii. 

As mentioned above, our simulation produces 
essentially two shocks located SE of the X-ray 
core whereas only one extended radio source has 
been observed that does not appear to be as- 
sociated with a particular galaxy. Both shocks 
appear to have a morphology, strength and evo- 
lution consistent with producing the SE radio 
source. That is, they are both smaller, more 
linear and slower moving than the NW shock 
while still being strong enough to produce sig- 
nificant particle acceleration. In terms of its lo- 
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cation with respect to the XSB, the SE shock 
furthest from the core is more consistent with 
the observed radio source. However, its prop- 
erties suffer from the same sensitivity to initial 
conditions as we described above regarding the 
NW shock. That said, the important point is 
that the single merger does produce shocks on ei- 
ther side of the cluster core which are capable of 
generating the observed radio morphology. Why 
some shocks within a cluster may produce ex- 
tended radio sources and others do not will likely 
depend on the specific details of the local mag- 
netic field structure and relativistic particle pop- 
ulations. These are the same issues that likely de- 
termine why some clusters have radio halos and 
others do not (see §|6l^). 

6.5. Why are there not more Halos of this 
type? 

If the large shocks are ubiquitous features of 
major merger events and mergers appear to be 
relatively common, why don't we see more ra- 
dio halos of this type? Specifically, according to 
our models, A754 is in a stage of evolution sim- 
ilar to (though somewhat earlier) than A3667, 
but it does not have the same type of radio 
halo (Roettiger et al. 1998a). Andernach et 
al. (1988) claims the existence of a radio halo in 
A754. However, higher resolution images indi- 
cate that this source is composed of three dis- 
crete sources, one of which is a NAT (Owen &; 
Ledlow 1997). The other two are likely back- 
ground objects. In any case, there does not ap- 
pear to be a halo in A754 that is comparable 
to A3667. There are several possible explana- 
tions for this apparent inconsistency. First, the 
shock structures, although long-lived relative to 
the canonical radio source age, are still relatively 
short-lived compared to the period between ma- 
jor mergers (2-4 Gyrs; Edge et al. 1992). Second, 
although the shocks may represent a necessary 
condition of halo formation, by themselves, they 
are not sufficient. Certainly, a large-scale mag- 
netic field is required and possibly a pre-existing 



population of relativistic electrons that are reac- 
celerated by the shock. At this time, very little 
is known about either the large-scale distribu- 
tion of magnetic field or relativistic particles as a 
general property of clusters. Some clusters may 
simply not have them or may have them in in- 
sufficient quantities. Third, cluster richness ap- 
pears to be a factor since only the richest appear 
to have halos (Hanisch 1982). Richer (i.e., more 
massive) clusters have deeper gravitational po- 
tentials which provide for more violent, higher 
velocity mergers and stronger shocks. They also 
have a more extended, higher pressure ICM, and 
more galaxies which may supply seed relativistic 
electrons and magnetic fields (Okoye & Onuora 
1997). Finally, projection effects may also play a 
role. 

7. SUMMARY 

We have presented a plausible model for the 
extended radio emission in A3667. In this model, 
A3667 has undergone a recent (~1 Gyr) merger 
event with a subcluster having a total mass of 
20% of the primary cluster. We believe that the 
subcluster impacted slightly off-axis while mov- 
ing from the SE to the NW. The merger gen- 
erates multiple shocks which provide sites for 
diffusive shock acceleration of relativistic elec- 
trons. The spatial extent of the radio sources 
is determined by the shock dimensions as well as 
the assumed strength and structure of the mag- 
netic field which ultimately determines the rate 
at which the relativistic particle age spectrally 
and the rate at which they diffuse/advect from 
the shock. This model is capable of reproduc- 
ing several features of the basic X-ray morphol- 
ogy including the twisted isophotes (although the 
twisting is more extreme in our model) , the steep 
X-ray surface brightness gradient to the SE (indi- 
cating a compression of the ICM), and the X-ray 
extension to the NW (indicating that the ICM 
has been drawn out of the cluster by the sub- 
cluster remnant). By modeling the spectral aging 
of the electron distribution, we estimate a mean 
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magnetic field of order less than 1 /xG. 

This model also explains many features of the 
A3667 radio morphology including: 1) the loca- 
tions of the NW and SE radio sources with re- 
spect to the X-ray surface brightness (these are 
the were the shocks are located) , 2) the sharp 
edge of the NW radio source (the site of particle 
acceleration), 3) the steepening of the spectral in- 
dex toward the southeast in the NW source (par- 
ticles age as they move away from the shock), 4) 
the relative sizes of the radio sources (the NW 
shock has a much greater lateral extent and is 
moving much faster than the SE shocks which are 
virtually stationary at this epoch), and, 5) the 
shapes of the radio components (the SE shocks 
are more linear than the NW shock). 

Although we have been able to reproduce the 
basic radio morphology, there are limitations to 
this model. The NW shock, in our model, is not 
strong enough to produce the flat (a = 0.5) spec- 
trum observed on the leading edge of the source. 
This shock could be strengthened if there is in- 
fall along an external filament (which we do not 
model) as suggested by Ensslin et al. (1998). 
The SE shocks are stronger than the NW shock 
largely because of residual infall from the sub- 
cluster which impinges on gas expelled from the 
cluster core by oscillations in the gravitational 
potential. It is also important to note that 
we do not evolve the relativistic particles self- 
consistently within the MHD simulation. Mod- 
eling the diffusion of particles from the shock is 
very complicated. For simplicity, we have as- 
sumed a weak field (high diffusion) limit in which 
particles are immediately left behind the shock. 
The true situation is undoubtedly more complex. 
Future work will entail evolving the relativistic 
particles in a more self-consistent manner as a 
second-fluid (e.g., Jun & Jones 1998; Jones et 
al. 1998). 

At this time we can only speculate as to why 
this type of radio source is so rare. It is most 
likely that very specific initial conditions are re- 
quired that may not exist in all clusters. It is 



possible that not all clusters contain a cluster- 
wide magnetic field or a seed population of rel- 
ativistic electrons. It is also likely that the rel- 
atively short radiative lifetime of the relativistic 
particles combined with the short lifetime of the 
shocks contributes to the rarity of these sources. 

We thank the Earth and Space Data Comput- 
ing Division at the NASA Goddard Space Flight 
Center for use of the Cray T3E supercomputer. 
We acknowledge support from NSF grant AST- 
9896039. We also thank H. Rottgering for use of 
Fig. 0. 

Note to be added in proof: Markevitch, 
Sarazin, & Vikhlinin (1998) recently produced 
a temperature map of A3667 using ASCA data. 
Our model is largely consistent with this map 
over the limited extent of the X-ray emission. 
The greatest inconsistency is the very hot gas 
SE of the X-ray core which is not apparent in 
the data. This would seem to indicate that our 
subcluster may be penetrating too deep into the 
cluster. 



12 



REFERENCES 

Andernach, H., Tie, H., Sievers, A., Reuter, 
H.-R, Junkes, N., & Wielebinski, R. 1988, 
A&AS, 73, 265 

Bagchi, J., Pislar, V., & Lima Neto, G. B. 

1998, MNRAS, submitted ( |astro-ph/9803020| ) 

Rinnpy. J. X^ TremRine, S. 1987, Oalartir Dy- 



659 

Fabian, A. C. 1994, ARAA, 32, 277 
Fadda, D., Girardi, M., Giuricin, G., Mardirossian, 
F., & Mezzetti, M. 1996, ApJ, 473, 670 

Feretti, L., Dallacasa, D., Giovannini, G., & 
Tagliani 1995, A&A 302, 680 



Feretti, L., &: Giovannini, G. 1997, astro- 



namicS) (Princeton: Princeton Universit^f 



ph/9709294 



Press) 
Biviano, A., Girardi, M., Giuricin, G., Mardirossian, 

F., & Mezzetti, M. 1993, ApJ, 411, L13 
Bryan, G. L., Klypin, A., Loken, C., Norman, 

M. L., & Burns, J. O. 1994, ApJ, 437, 5L 

Burns, J. O., Roettiger, K., Pinkney, J., Per- 

ley, R. A., Owen, F. N., & Voges, W. 

1995, ApJ, 446, 583 
Burns, J. O., Roettiger, K., Ledlow, M. & 

Klypin, A. 1994, ApJ, 427, L87 
Burns, J. O. et al. 1997, in Galactic & Cluster 

Cooling Flows, ed. N. Soker (ASP:San 

Francisco), Vol. 115, 21 
Burns, J. O. 1998, Science, 280, 345 
Crone, M., Evrard, A., k Richstone, D. 1994, 

ApJ, 434, 402 
DeYoung, D. S. 1992, ApJ, 386, 464 

Deiss, B. M., Reich, W., Lesch, H., & Wielebin- 

ski, R. 1997, A&A, 321, 55 
Dressier, A. & Shectman, S. A. 1988, AJ, 95, 

985 
Drury, L. 1983, Rep. Prog. Phys., 154, 973 
Edge, A. C, Stewart, G. C., & Fabian, A. C. 

1992, MNRAS, 258, 177 
Eilek, J. A. & Henriksen, R. N. 1984, ApJ, 

277, 820 
Ensslin, T. A., & Biermann, P. L. 1998, A&A, 

330, 90 
Ensslin, T. A., Biermann, P. L., Klein, U., & 

Kohle, S. 1998, A&A, 332, 395 

vrarA A F 1 QQn Ap.T .'^fi.'^ ."^^Q 



I 



Goldman, I. & Rephaeli, Y. 1991, ApJ, 380, 
344 

Gomez, P. L, Loken, C, Roettiger, K. & 
Burns, J. O. 1999, in preparation 

Goss, W. M., Ekers, R., Kellern, D. J., k 
Smith, R. M., 1982, MNRAS, 198, 259 

Gropp, W., Lusk, E., k Skjellum, A. 1994, 
Using MPI: Portable Parallel Program- 
ing with the Message-Passing Interface, 
(Cambridge: MIT Press) 

Hanisch, R. J. 1980, AJ 85, 1565 

Hanisch, R. J. 1982, A&A 116, 137 

Hawley, J. F., k Stone, J. M. 1995, Comp. 
Phys. Comm., 89, 127 

Hockney, R., k Eastwood, J. 1988, Computer 
Simulation Using Particles, (Philadelphia: 
lOP) 

Holman, G. D., lonson, J. A., Scott, J. S. 1979, 
ApJ, 228, 576 

Hunstead, R. W. k Wieringa, M. 1998, private 
communication. 

Jackson, J. D. 1975, Classical Electrodynam- 
ics, (New York: Wiley) 

Jaffe, W. J. 1977, ApJ, 212, 1 

Jaffe, W. J. k Perola, G. C. 1973, A&A, 26 
423 

Jones, P. A., k McAdam, W. B., 1992, ApJS, 
80, 137 

Jones, T. W., Ryu, D., k Engel, A. 1998, ApJ, 



astro-ph/9809081 



Jun, B.-I. k Jones, T. W. 1998, ApJ, astro- 



vans, C. R. k IIawle)>, J. F. 1988, ApJ, 332, 



ph/9809082 



Kardashev, N. S. 1962, Soviet Astr. -AJ, 6, 



13 



317 
Kennel, C. F., Coroniti, F. V., Scarf, F. L., 

Livesey, W. A., Russell, C. T., Wenzel, E. 

J. & Scholer, M. 1986, J. Geophys. Res. 

91, 11 
Kim, K.-T., Kronberg P. P., Dewdney, P. E., 

Landecker, T. L., 1990, ApJ, 355, 29 
Kim, K.-T., Kronberg, P. P., Tribble, P. 1991, 

ApJ, 379, 80 

Knopp, G. P., Henry, J. P., & Briel, U. G. 

1996, ApJ, 472, 125 
Longair, M. S. High Energy Astrophyics, (Cam' 

bridge, UK: Cambridge Univ. Press) 
Markevitch, M., Forman, W. R., Sarazin, C. 

L., k Vikhlinin, A. 1998, ApJ, 503, 77 

Markevitch, M., Sarazin, C. L., & Vikhlinin, 
A. 1998, ApJ Letters, submitted ( [astro- 
ph/9812005| ) 



McGlynn, T. A. & Fabian, A. C, 1984, ApJ, 

208, 709 
Miley, G, 1980, ARAA, 18, 165 
Mohr, J. J., Evrard, A. E., Fabricant, D. G., 

& Geller, M. J. 1995, ApJ, 447, 8 

Myers, S. T., & Spangler, S. R. 1985, ApJ, 

291, 52 
Navarro, J. F., Frenk, C. S., & White, S. D. 

1995, MNRAS, 275, 720 
Navarro, J. F., Frenk, C. S., & White, S. D. 

1997, ApJ, 490, 493 
Okoye, S. E. & Onoura, L. I. 1997, MNRAS, 

283, 1047 
Owen, F. N,. & Ledlow, M. J. 1997, ApJS, 

108, 41 

Pacholczyk, A. G. 1970, Radio Astrophysics, 
(San Francisco: Freeman) 

Rottgering, H., Snellen, I., Miley, G., de Jong, 

J. P., Hanisch, R., & Perley, R. A. 1994, 

ApJ, 436, 654 
Rottgering, H., Wieringa, M. H., Hunstead, 

R. W., & Ekers, R. D. 1997, MNRAS, 

290, 577 



Roettiger, K., Burns, J. O., & Loken, C 1993, 

ApJ, 407, 53L 
Roettiger, K., Burns, J. O., & Loken, C 1996, 

ApJ, 473, 651 
Roettiger, K., Loken, C, & Burns, J. O. 1997, 

ApJS, 109, 307 
Roettiger, K., Stone, J. M., & Mushotzky, R. 

F. 1997, ApJ, 482, 588 
Roettiger, K., Stone, J. M., & Mushotzky, R. 

F. 1998a, ApJ, 493, 62 
Roettiger, K., Stone, J. M., & Burns, J. O. 

1998b, ApJ, in press. 

Roland, J. 1981, A&A, 93, 407 
Ruzmaikin, A., Sokoloff A., & Shukurov, A. 
1989, MNRAS, 241, 1 

Sarazin, C. L. 1988 in X-ray Emissions from 
Clusters of Galaxies. Cambridge Univer- 
sity Press, Cambridge 

Schilizzi R. T., & McAdam, W. B 1975, Mem. 
R. Astro. Soc, 79, 1 

Sodre, L., Capelato, H. V., Steiner, J. E., 
Proust, D., Mazure, A. 1992, MNRAS, 
259, 233 

Stone, J. M. & Norman, M. 1992a, ApJS, 80, 
753 

Stone, J. M. & Norman, M. 1992b, ApJS, 80, 
791 

Stone, J. M., Hawley, J, F., Evans, C. R., & 
Norman, M. L. 1992, ApJ, 388, 415 

Tao, L. 1995, MNRAS, 275, 965 

Tribble, P. C. 1991 MNRAS, 253, 147 

Tribble, P. C. 1993a MNRAS, 261, 57 

Tribble, P. C. 1993b MNRAS, 263, 31 

Tribble, P. C. 1994 MNRAS, 269, 110 



14 



FIGURE CAPTIONS 

Fig. I The ROSAT FSPC X-ray image (0.1- 
2.4 keV, contours) overlaid on the grayscale of 
the wide-field, 843 MHz MOST image. The X- 
ray contours correspond to 2, 8, 18, 32, 50, 72, 
98, 128, 162, 200, and 242 times the background 
noise. The rms noise level of the radio image is 
0.7 mJy beam^^. (Figure courtesy H. Rottgering 
from Rottgering et al. 1997). 

Fig. pi Galaxy positions (both filled and 
open circles) overlaid on the ROSAT PSPC im- 
age (contours). The circle sizes are proportional 
to the galaxy's 6 statistic (Dressier & Shectman 
1988). A large value of 6 indicates a large local 
deviation in velocity mean and dispersion from 
the global values. We have filled circles repre- 
senting local distributions having 5 less than 20% 
of the maximum. Note the two regions associated 
with the XSB peak and to the NW of the peak. 
Both of these are associated with D galaxies. 

Fig. 1^ a) The simulated X-ray surface bright- 
ness (dashed contours) and projected emission- 
weighted temperature (solid contours). We have 
added noise to the synthetic X-ray image which 
has been contoured at the same levels (relative 
to peak) as the ROSAT image (see Fig. 0). The 
temperature contours are labeled with the corre- 
sponding temperature in keV. The image is 3.15 
X 3.85 Mpc. b) The ASCA temperature profile ( 
Markevitch et al. 1998, open diamonds) and the 
simulated temperature profile (asterisks) within 
the same spatial bins (as indicated by the hor- 
izontal bar, assuming Ho=70 km s^^ Mpc^^). 
The profiles have been slightly displaced for clar- 
ity. The vertical bars on the ASCA data indicate 
the uncertainty in the temperature while the ver- 
tical bars on the simulated temperatures repre- 
sent the dispersion in temperature withing the 
volume. Note that although the cluster is non- 
isothermal, it appears to be quite isothermal in 
profile. 

Fig. ^ a) The simulated A3667 X-ray sur- 
face brightness and radio data. Contours rep- 



resent the X-ray surface brightness. We have 
added noise to the synthetic X-ray image which 
has been contoured at the same levels (relative 
to peak) as the -ROS'^T image (see Fig. ||). The 
grayscale represents the synthetic radio emission 
at 1.4 GHz. The + symbols indicate the loca- 
tions of peaks in the dark matter /galaxy distri- 
bution. The dashed and dotted lines traversing 
the synthetic radio emission to the NW (upper 
right) indicate the location of the radio spectral 
index {a\-Q) profiles displayed in (b). The dashed 
and dotted lines in (a) correspond to the dashed 
and dotted lines in (b). The profiles extend from 
the southeast (lower left) to the northwest (upper 
right). The image is 3.15 x 3.85 Mpc. 

Fig. |5| The simulated X-ray surface bright- 
ness (contours, same as in Figs. ^ and ^) over- 
laid onto the underlying shock structure (solid 
shading) within the plane of the merger. Regions 
with a density compression greater than two are 
shaded black. The image is 3.15 x 3.85 Mpc. 
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Table 1. Initial Cluster Parameters 



Cluster^ 
ID 


ML 
(1014 Mo) 


(keV) 


(km s-i) 


(kpc) 


Jg 


v' 

"impact 

(km s-i) 


1 

2 


13.0 
2.6 


6.5 
3.6 


943 
540 


288 
144 


0.045 
0.045 


2750 



1 Cluster ID. ^ Total Mass R<3 Mpc. ^ Temperature. ^ 1-Dimensional Velocity Dispersion R<1.5 Mpc. 
^ Core Radius. ^ Global Gas Fraction, by mass. ^ Impact Velocity (dark matter). 
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